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Abstract 

q ■ We discuss the spatiotemporal intermittency (STI) seen in coupled map 
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lattices (CML-s). We identify the types of intermittency seen in such sys- 
tems in the context of several specific CML-s. The Chate-Manneville CML 
is introduced and the on-going debate on the connection of the spatiotempo- 



' ral intermittency seen in this model with the problem of directed percolation 

(N ' 

is summarised. We also discuss the STI seen in the sine circle map model 

o 
in 

■ mogenous logistic map lattice which shows the novel phenomenon of spatial 



and its connection with the directed percolation problem, as well as the inho- 



intermittency and other types of behaviour not seen in the other models. The 
connection of the bifurcation behaviour in this model with STI is touched 



•l-H 

15 

> 

^ , upon. We conclude with a discussion of open problems. 

1 Introduction 

Regimes of spatiotemporal intermittency, wherein laminar regions, which exhibit 
regular dynamics in space and time, coexist and propagate together with regions 
which show irregular or chaotic bursts, can be seen in a wide variety of spatially 
extended systems in the laboratory as well as in theoretical model systems. Such 
regimes have been regarded as precursors to regimes of fully developed spatiotem- 
poral chaos p. Spatiotemporal intermittent behaviour has been seen in theoretical 
models such as coupled map lattices pOE], probabilistic cellular automata [I], and 



partial differential equations |5J E] , as well as in experimental systems such as chem- 
ical reactions [7j, Rayleigh-Benard convection in narrow channels and annuli [HUH]; 
planar Couette flow [TU], studies of fluid flows between rotating eccentric cylinders 
such as the Taylor-Dean ^Tj and Taylor- Couette |T21 ITS] flows, and the 'printer's 
instability' 14 . A variety of scaling laws have been observed in these systems. 
However, there are no definite conclusions about their universal behaviour. Many of 
the observed phenomena have been seen in experimental systems where no simple 
model is available. There is a lively on-going debate about the nature of spatiotem- 
poral intermittency and its analogy with systems which undergo phase transitions. 
Thus, the analysis of spatiotemporal intermittency remains a challenging theoreti- 
cal problem. In this review, we will concentrate on spatiotemporal intermittency in 
coupled map lattices, i.e. systems with continuous variables which evolve on discrete 
space-time. 

Spatiotemporally intermittent behavior has been divided into two classes. In the 
first type, there is no spontaneous creation of bursts. Under spatial coarse-graining, 
if a site and its neighbours show laminar behaviour at a given step, the site remains 
laminar at the next time step. However it can show turbulent behaviour if infected 
by a turbulent neighbour. A turbulent site can either relax to laminar behaviour or 
infect its neighbour. Below a threshold parameter, which controls the propagation 
rate of turbulence, the turbulence dies out. Once the whole system reaches the 
laminar state, it has no way of escaping from this state. Hence the laminar state is 
called the absorbing state. Above the threshold parameter, the turbulence spreads 
over the whole lattice. Thus the burst states 'percolate' through space-time. It has 
therefore been argued that the transition to this type of spatiotemporal intermit- 
tency is a second order phase transition, and that this transition falls in the same 
universality class as directed percolation [T^j with the laminar states being identified 
with the 'inactive' states and the turbulent states being identified as the 'active' or 



percolating states. This conjecture has become the central issue in a long-standing 
debate [TBJ 113 CHI EE [2D] , which is still not completely resolved. 

In the second type of intermittency, the spontaneous creation of bursts can be 
seen, given some coarse-grained reduction of the states. Thus, for this laminar 
site has a finite probability of undergoing a turbulent burst even if all its neighbours 
are laminar. This kind of intermittency is associated with the transition of a spatial 
pattern. A typical example of this type of intermittency is observed in the coupled 
logistic lattice from pattern selection to fully developed spatiotemporal chaos j2j as 
well as in experiments on Benard convection [HIE], Faraday instabilities [21j . and in 
two-dimensional electric convection in a liquid crystal [22] • 



2 The C hate- Marine ville map and its variants: 



The conjecture by Pomeau that the transition to spatiotemporal intermittency of 
the first type with a laminar absorbing state belongs to the same universality class as 
directed percolation has been very extensively investigated. A very simple coupled 
map model was introduced by Chate and Manneville for the investigation of this 
conjecture. 

The Chate- Manneville maps at different sites are coupled together by linear 
diffusive coupling of the form 



x «+i(i) = (i _ e)f(x n (i)) + e -[f{x n {i + 1)) + /(*»(* - 1))] 



(1) 



The single site Chate- Manneville map has the form 
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where r > 2. The return map is shown in Fig. 1. 

The uncoupled dynamics of such a map is chaotic as long as / (x) remains in the 
unit interval since r > 1 everywhere in this domain. However, as soon as f(x) > 1 
the iteration is locked on to a fixed point. The local phase space is thus a marginally 
stable fixed line ( x G [l,r/2]), connected to a chaotic repeller ( x G [0, 1]). There 
is thus a clear and natural distinction between laminar and turbulent states for a 
given site. Chate and Manneville argued that since the model has no single stable 
fixed point or globally attracting period, there is no additional local dynamics to 
complicate the transition from local transient chaos to global spatiotemporal chaos 
via the spatial coupling alone. 

The dynamics of this chain was studied with random initial conditions. At 
e > e c , sustained spatiotemporal intermittency was seen for the system for r = 2.1 
and r = 3. The lattice showed regions of laminar behaviour punctuated by chaotic 
bursts (See Fig. 2). The distribution of the lengths of the laminar regions showed 
power-law behaviour of the form P(l) « with the exponent £ taking values 
between 1.67 and 1.78 at r = 2.1 and 1.9 to 1.99 at r = 3. (The simulated DP values 
for the corresponding exponent ranged from 1.6 to 1.74). It was initially argued 
that the values obtained showed significant departures from the directed percolation 
exponents and hence that spatiotemporal intermittency did not belong to the DP 
universality class. However Grassberger et al examined the same CML later 19\ and 
found that the departures from DP exponents at r = 2.1 were insignificant on long 
time-scales. In contrast, statistically significant departures from DP were observed 
at r = 3.0. It was observed that for this case, long straight lines of turbulent sites 
were seen in the space-time plot (Fig. 2(b)). These were called 'solitons' and were 
reflective of long range correlations in the system. Grassberger and Schreiber argued 
that these long range correlations spoilt the analogy with DP and hence led to the 
departure of the STI in this system from the DP universality class. The existence 



of solitons led to the emergence of a length and time scale (the mean life time of the 
solitons) which is independent of critical fluctuations, and hence led to departures 
from the DP class JH]- Since the solitons in this model had very long life-times, they 
also conjectured that the cross-over times in this model could be very long and the 
real asymptotic behaviour of the model could still be in the DP universality class. 

The conjecture that the solitons were the primary ingredient which spoilt DP 
universality was tested in a variety of ways [To] 123] 121] . The CML of eq. |2] was 
evolved with asynchronous up-dates i.e. the maps at every site were not updated 
simultaneously, but were updated one after another in a random sequence. The 
asynchronous updates served to destroy the solitons and good agreement with a 
complete set of DP exponents was obtained (see Appendix for definitions of the DP 
exponents and Table 1 for obtained values). 

Another study introduced a deterministic extension of the Chate-Manneville 
model, where the soliton properties could be tuned [2H1 El] • It was seen that the 
role played by the solitons in the transition to turbulence was more profound than 
that of setting a cross-over time. They were able to change nature of the transition 
from second to first order. 

The CML equations of this solitonic Chate-Manneville model were 

x n+1 (i) = (l-e)f(x n (i)) + ~[f(x n (i + l)) + f(x n (i-l))]+y n (i) 
y n+1 (i) = b{x n+1 {i) -x n (i)) 

(3) 

where the function f{x) is the same as that in Eq. El The map becomes increasingly 
two dimensional with increase in |6|. Thus if x* is a fixed point of the old CML Eq. 
121 it gets mapped to the fixed point (x*,0) of Eq. |3] The CML thus has three 
adjustable parameters r, e, b. The parameter r was set to 3.0 and the phase diagram 
in e, b space was obtained. The parameter b has a large influence on the existence of 



solitons. The 6 = case is the original Chate-Manneville case where there are long 
lived solitons (Fig. 3 (c),(d)). These life-times become even longer when b = —0.1 
(Fig. 3(a),(b)), and the solitons only vanish when they collide with other solitons or 
propagate into turbulent structures (Fig. 3(e), (f)). The order parameter m which is 
the number of active or turbulent sites, jumps sharply as a function of e — e c as in a 
first order transition [231 [21! • The behaviour of m at b = 0.0 and b = 0.2 is consistent 
with a continuous transition with exponents that approach DP values more closely 
as \b\ increases corresponding to soliton-poor regimes (Table. 2). Thus not only does 
the existence of solitons spoil the analogy of spatiotemporal intermittency with DP, 
but it can change the order of the transition itself. Thus, in order to obtain a 
clean comparison between directed percolation and spatiotemporal intermittency, it 
is desirable to examine a deterministic system which is free of coherent propagating 
structures. We describe such a system in the next section. 



3 The sine circle coupled map lattice: 

The sine circle map lattice is defined by 

6 n+1 (i) = (l-e)f(6 n (i)) + ~[f(e n (i + l))+f(e n (i-l))] (mod 1) (4) 

with periodic boundary conditions. The local map is the sine circle map 

f(9 n ) =e n + tt-— sin(27T0 n ) (5) 

Here, K is the strength of nonlinearity, Q is the frequency of each single sine 
circle map and, e is the coupling strength. This system is a popular model for the 
behaviour of mode-locked oscillators. The system shows regimes of spatiotemporal 
intermittency at certain points in the parameter space when random initial con- 
ditions are evolved synchronously. The phase diagram of the system at K = 1 is 



shown in Fig. 4 and the points where spatiotemporal intermittency are seen are 
marked with diamonds [20] and pluses [23]. This system shows a transition to spa- 
tiotemporal intermittency in the presence of two types of absorbing states. The two 
qualitatively distinct absorbing regions are as follows [20] . 

(i) When the nonlinearity parameter K = 1, there are regions of (e — Q) space 
where the system goes to the synchronised spatiotemporal fixed point 6* = ^- sin^^p). 
This constitutes an unique absorbing state. The critical behaviour was examined at 

2 critical points in this regime: Q = 0.064, e = 0.63775, and Q = 0.068, e = 0.73277. 
These mark the transition from a laminar phase to STL The turbulent sites here 
are those which are different from 9*. 

(ii) When the nonlinearity parameter K = 3.1, there are regions of (e — Q) space 
where sites with any value less than 1/2 constitute the absorbing states, and sites 
whose values are greater than 1/2 constitute the turbulent states. The absorbing 
states are now infinitely many, as also weakly chaotic. The critical behaviour was 
examined at 2 critical points in this regime as well: Q = 0.18, e = 0.701, and 
tt = 0.19, e = 0.65612. 

The exponents for the onset of spatiotemporal intermittency have been defined 
both for the bulk critical exponents and the spreading exponents (See Tables 3 and 
4). In both cases there is good agreement with DP values [22 . Thus the coupled 
sine circle map constitutes a deterministic system evolving under synchronous up- 
dates where DP exponents are straight-forwardly obtained. There are no long-lived 
coherent structures at any of the parameter values where the system has been studied 
to spoil the analogy with DP. This also lends credence to the argument that it is the 
solitons which lead to the departure of the Chate-Manneville exponents from the 
DP class. Certain parameter regions of the sine circle map model show regimes of 
spatiotemporal intermittency which is unlike the DP behaviour. These are currently 
under investigation. Such regimes have been seen earlier in a coupled map model 



which we will discuss below. 



4 The inhomogeneous logistic map lattice: 

The specific system studied by us consists of an inhomogeneous lattice of coupled 
maps where the inhomogeneity appears in the form of different values of map pa- 
rameter at different sites. Such lattices have been considered in the case of pinning 
studies j2H]- The coupling between maps is diffusive and nearest neighbour. We 
consider the simplest case, that is, a situation where neighbouring lattice sites have 
different values of the parameter and alternate sites have the same value of the 
parameter. 

This model is defined by the evolution equations 

x n+1 (2i) = (1 - e)f(x n (2t}^) + ~[f(x n (2i - 1),//) + f(x n (2i + 1),//)] (6) 

where f(x) = fix(l — x) is the logistic map and /z, /i' G [0,4], x n {2i) is the value 
of the variable x at the even lattice site 2i at time n, and < x < 1. In the case of 
x n (2i + 1), the variable defined at odd lattice sites, the evolution equation is defined 
by the evolution above with 2i being replaced by 2i + 1 and /i and interchanged. 
We set y! = fi — 7 and use periodic boundary conditions where 2N, the number of 
lattice sites is even. The synchronised fixed points of the system are given by x* = 
and x* = H=2L=1. 

Spatiotemporal intermittency can be found at many points in the parameter 
space, unlike the homogeneous logistic coupled map lattice, where STI is found for 
values of fi which lie near the period 3 window. 

Bifurcations from the synchronised fixed points can be used to obtain indica- 
tors of regions where spatiotemporal intermittency can be found. Spatiotemporal 
intermittency can be seen in the vicinity of bifurcation lines of co-dimension 1 and 



bifurcation points of co-dimension 2 where such lines intersect each other. It is in- 
teresting to note that pure spatial intermittency, where there are bursts interspersed 
with laminar behaviour on the spatial lattice, accompanied by purely temporal peri- 
odic behaviour, can be found in the vicinity of some tangent-period doubling points 
where the eigen- values of the stability matrix cross +1 and —1 at the same point in 
the parameter space. 

These conditions are satisfied for the fixed point x* = 0.0 in the neighbourhood 
of e = 0.635, 7 = 1.1 where the eigen- values are —0.9127, and 1.0585 respectively. 
We show the spatial intermittency present in the vicinity of this point in Fig. 5. 
The parameter values chosen are 7 = 1.19 and e = 0.635. The space-time plot of 
the system can be seen in Fig. 5. The temporal period two structure of both the 
laminar and burst sites can be seen clearly in the space-time plot. It is clear that 
the time evolution of the lattice shows stable period two oscillations at both the 
laminar and the burst sites. Thus the intermittency is purely spatial in nature and 
temporally we have stable periodic behaviour. 

The temporally periodic spatial structure seen here shows the unusual occurrence 
of long range spatial correlations. This is reflected in the distribution of laminar 
lengths. The length of the laminar lengths, i.e. the number of consecutive sites 
which follow periodic behaviour before being interrupted by bursts is calculated. As 
mentioned above the spatial period in the laminar regions is two. The distribution 
of laminar lengths shows power-law behaviour P(l) m with exponent ( = 1.1. 

Several other tangent-period-doubling points are also seen in this parameter 
space and spatial intermittency with an associated scaling exponent ( = 1.1 can 
be seen in their vicinity , except for the tangent-period doubling point at 7 = 0.66 
and e = 0.39 where spatiotemporal intermittency is seen (See Fig. 6(a), (b)). The 
distribution of laminar lengths seen here also shows power-law scaling P(l) = l~^ F 
but now with exponent (p = 2.0. Spatiotemporal intermittency can also be seen 



along lines of co-dimension 1 bifurcations along the parameter space. Here, other 
exponents (2 = 1-33 and £3 = 1.66 are seen. These appear to be crossover exponents. 

It is interesting to note that an exponent which falls within the same range as (f 
has been seen in the case of Rayleigh-Benard convection in an annulus exponents 
that fall in the same range as (2 have been seen for convection in a channel and for 
the Taylor-Dean experiment, and £3 has been seen in the Chate-Manneville CML 
at r = 2.1 and e = 0.0045 as well as in DP simulations ^2]- Similarly the Chate- 
Manneville CML at r = 3.0 e c = 0.360 shows the exponent (p. The behaviour in 
this model is complicated by the fact that the laminar state is periodic, and that 
a bifurcation to period two seems to take place very close to the bifurcation from 
the synchronised fixed point, therefore it is difficult to draw direct connections to 
behaviour seen in the earlier models. However, the model provides pointers to the 
following conjectures. Spatiotemporal intermittency can be found in the vicinity of 
bifurcations from low-ordered periodic points. The type of bifurcation that takes 
place as well as the nature of the stable state which precedes the bifurcation may 
provide an indication of the universality class of the spatiotemporal intermittency. 
Further work on this system could lead to interesting results. 

5 Discussion 

To summarise, spatiotemporal intermittency is a dynamical regime with interesting 
but little understood behaviour. While a plethora of phenomenological information 
is available on the behaviour of systems which show STI, both theoretical and ex- 
perimental, a unifying picture is lacking. While the conjecture that several systems 
which demonstrate spatiotemporal intermittency belong to the DP universality class 
finds some support, further work is necessary to identify the factors which contribute 
to the presence or absence of DP. 



The effect of travelling coherent structures on the transition to spatiotemporal 
intermittency deserves to be explored further, as does the connection with bifur- 
cation behaviour. Studies of systems where the laminar state is not absorbing are 
few and far between. Systematic studies of such systems need to be undertaken. 
Further work on experimental systems is also necessary. There has been only one 
experiment where DP exponents have been found so far [3T]. Further work also 
needs to be done in constructing models for the quasi- Id fluid systems which show 
spatiotemporal intermittency. The extent to which the behaviour seen in stochastic 
models like probabilistic cellular automata, is similar to the behaviour in strictly 
deterministic models, like those considered here, also needs to be explored further. 
Thus, we expect this area of research to remain active and contribute to many lively 
debates in the coming years. 
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6 Appendix 

The quantities characterised by critical exponents at e = e c , the critical value at 
which the transition to sustained spatiotemporal intermittency takes place, are : 

1. The static exponents: 

(a) The escape time, r, i.e. the time taken by the system, starting from 
a random initial state, to reach the absorbing state, shows power law 
behaviour of the form r ~ L z . 

(b) The order parameter m, i.e. the fraction of turbulent sites, scales as 
m(e c , L, t) ~ (e — e^ 13 , for e > e c , at fixed t. 

(c) For t <^ r, m (e c ,L,t) ~ t-M" z , where v is the correlation length exponent 
in space. 

(d) The correlation function Cj(t) is defined as 

1 L 

C M = J E < Xi(t)x i+j (t) >-< Xi (t) > 2 (7) 
L i=i 

where Xi(t) is the variable value at site % at time t. The brackets denote 
averaging over different initial conditions. This quantity scales as Cj (t) ~ 

(e) These bulk exponents satisfy the hyperscaling relation 2(3 /v = 2 — d + rj'. 

2. When the turbulent seeds spread in an absorbing lattice, the quantities asso- 
ciated with the spreading exponents are : 

(a) N(t), the number of active sites at time t. This scales as N(t) ~ t v , 

(b) P(t), the survival probability, or the fraction of initial conditions which 
show a non-zero number of active sites at time t, which shows scaling 
behaviour of the form P(t) ~ t~ 5 . 



(c) R 2 (t), the mean squared deviation of the position of the active sites from 
the original sites of turbulent activity, averaged over the surviving runs 
alone. This scales as R 2 {t) ~ t Zs 
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Table 1: The critical exponents of the asynchronously updated Chate-Manneville 
coupled map lattice. The last row shows the corresponding DP exponents. 
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Table 2: The critical exponents of the extended 2 dimensional Chate Manneville 
CML for different values of b. b = 0.2 shows exponents closer to the DP exponents. 
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Table 3: Critical static exponents of the synchronously updated coupled circle map 
lattice for 4 critical points. The first two critical points correspond to transitions to 
an unique absorbing case, while the third and fourth points correspond to weakly 
chaotic absorbing states. The last row shows the corresponding exponents of directed 
percolation. 
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Table 4: Spreading exponents of the synchronously updated coupled circle map 
lattice for 4 critical points. Two active seeds in an absorbing configuration is used 
as initial condition. For the first 2 critical points there exists an unique absorbing 
state, while for the third and fourth points one can have many different absorbing 
states and consequently many different initial absorbing backgrounds. However we 
notice that the exponents obtained are the same for different initial preparations 
and thus appear universal for 2 active seeds. The last row shows the corresponding 
exponents of directed percolation. 




Figure 1: The Chate-Manneville map at r = 2.5 
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(a) 




0.0045 and (b)r = 3.0, e c = 0.36 after discarding 3200 transients. Laminar sites are 
the black spaces. In (a), every 8th timestep has been plotted. The slower dynamics 
at r = 2.1 can be seen from the plot. Also, the solitons are clearly seen in (b). 




Figure 3: Space-time plots at points above (left column) and below (right column) 
criticality for the 2 dimensional Chate-Manneville map for different values of b. (a) 
b = -0.1, e = 0.359 , (b) b = -0.1, e = 0.3537, (c) b = 0.0, e = 0.370 , (d) b = 0.0, 
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Figure 4: The phase diagram in (e, Q) space for the coupled circle map lattice at 
K — 1 after discarding 15000 transients for a lattice of 1000 sites. Random initial 
conditions were used. The regimes where spatially synchronised and temporally fixed 
solutions are obtained are shown by dots. The points at which DP-like exponents 
are obtained are marked by O [20] and + [25]. 
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Figure 5: Spacetime plot of the inhomogenously coupled logistic map at 7 = 1.19 
and e = 0.635 for a system of 500 sites. The last 25 iterations are shown after 
discarding 5000 transients. Pure spatial intermittency, where the laminar region is 
temporally periodic is seen. 




Figure 6: Spatiotemporally intermittent behaviour seen in the inhomogenously cou- 
pled logistic map at 7 = 0.66 and e = 0.39. (a) L = 2000, last 100 iterations are 
shown after discarding 20, 000 transients, (b) L = 500. The last 25 iterations are 
shown after discarding 5000 transients. 



